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ABSTRACT 


The  U.S.  Department  of  Energy's  (DOE)  procedure  for  testing  the  Annual  Fuel  Utilization 
Efficiency  (AFUE)  of  residential  central  furnaces  and  boilers  references  ASHRAE  Standard  103- 
1993.  The  standard  has  been  used  by  the  fumace/boiler  industry  to  obtain  the  AFUE  values  of 
their  products  as  required  by  DOE.  ASHRAE  Standard  Project  Committee  (SPC)  103R  is  in  the 
process  of  making  revisions  to  the  standard  to  address  some  shortcomings  and  to  improve  the 
calculation  procedure  on  the  basis  of  comments  and  experiences  from  the  industry.  This  report 
presents  proposed  revisions  to  certain  part  of  the  test  procedure  to  address  some  of  the 
comments.  The  proposed  revisions  include  a new  procedure  to  calculate  the  average  burner 
on/off  times  of  modulating  furnaces/boilers  during  the  cyclic  operation  of  the  burner,  a revised 
method  for  the  calculation  of  the  annual  burner  operating  hours,  and  an  extension  of  the 
procedure  to  calculate  the  off-period  flue  loss  of  oil-fired  equipment  designed  with  a long  post 
purge  time.  This  report  also  provides  the  background  information  on  the  derivation  of  the 
integral  functions  listed  in  Appendix  E of  ASHRAE  Standard  103-1993  that  have  not  been 
formally  documented.  Those  integral  functions  were  used  for  producing  the  graphs  in  the 
standard  for  the  hand  calculation  of  the  off-period  losses  with  a calculator,  and  also  in  the 
numerical  algorithm  in  the  computerized  calculation  procedure. 
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steady-state  efficiency  in  percent 

average  steady-state  efficiency  operating  in  modulating  mode  in  percent 
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higher  heating  value  of  fuel 
heating  load  hours 

multiplication  factor  for  sensible  heat  loss  during  the  off-period 
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off-cycle  sensible  heat  loss  in  percent 

off-cycle  flue  gas  mass  flow  rate  (per  minute) 
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electrical  power  (kW)  to  burner 

steady-state  rated  (maximum)  fuel  input  rate  (per  hour) 

steady-state  maximum  fuel  input  rate  (per  hour)  for  a modulating  unit 

steady-state  reduced  (minimum)  fuel  input  rate  (per  hour) 

heating  capacity,  useful  heat  output  per  hour,  under  steady-state  rated  fuel  input 

average  heat  output  rate  (per  hour)  for  modulating  mode 

heating  capacity,  at  reduced  (minimum)  fuel  input  rate 

fuel  energy  input  per  hour  to  pilot  light 

time  in  minute 

power  burner  post-purge  time  period 

burner  off-time  per  cycle 

burner  on-time  per  cycle 

balance  point  temperature  for  modulating  units 

design  outdoor  air  temperature 

off-period  flue  gas  temperature  at  time  t 

steady-state  flue  gas  temperature 
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1.  INTRODUCTION 


The  U.S.  Department  of  Energy  (DOE)  test  procedure  [l]  for  determining  the  Annual  Fuel 
Utilization  Efficiency  (AFUE)  of  residential  central  furnaces  and  boilers  referenced  the  industry 
test  standard  ANSI/ASHRAE  103-1993  [2]  (ASHRAE  103-1993)  of  the  American  Society  of 
Heating,  Refrigerating,  and  Air  Conditioning  Engineers  (ASHRAE).  The  standard  was  initially 
developed  in  1982  based  on  the  DOE  test  procedures  for  single  stage  furnaces  and  boilers 
recommended  by  Kelly  et  al [3].  It  was  subsequently  revised  in  1988  and  again  in  1993  to  include 
test  procedures  for  condensing  units,  modulating  (or  multi-fire)  units,  and  for  units  employing  a 
short  post-purge  period  after  the  burner  is  shut  off.  In  1998  ASHRAE  organized  the  Standard 
Project  Committee  (SPC)  103R  to  begin  the  revision  process  to  ASHRAE  103-1993  following 
comments  from  the  industry  on  the  need  to  address  some  possible  shortcomings  of  the  standard 
based  on  user  experiences.  This  report  documents  the  proposals  to  SPC-103R  to  address  the 
issues  raised  by  industry  concerning  the  test  standard.  The  issues  are  as  following: 

(1)  The  inconsistency  in  the  calculated  AFUE  values  for  a two-stage  furnace  (or  boiler) 
when  they  are  calculated  on  the  basis  of  (a)  the  single-stage  procedure  and  (b)  the  modulating 
two-stage  procedure  using  the  same  test  data.  The  first  one  treats  the  furnace  as  a single  stage 
furnace  operating  at  the  reduced  input  to  supply  the  heating  load,  and  the  second  one  treats  the 
furnace  as  two-stage  furnace.  With  a ratio  of  reduced  input  to  rated  input  of  from  60  % to  70  %, 
the  standard  shows  that  the  furnace  would  supply  almost  98  % of  the  annual  heating  load  at  the 
reduced  input.  That  is,  the  furnace  is  operating  as  a single-stage  furnace  at  the  reduced  fuel  input 
rate  most  of  the  time  even  when  it  is  treated  as  a two-stage  furnace.  However,  the  calculated 
results  showed  a difference  in  the  AFUE  values  of  more  than  one  percentage  point  between  the 
two  procedures. 

(2)  In  the  calculation  for  the  number  of  annual  Burner  Operating  Hours  (BOH)  for  a 
modulating  furnace  (or  boiler),  the  standard  estimates  the  BOH  using  the  fuel  and  auxiliary 
electric  energy  input  measured  at  the  maximum  rated  input  condition.  For  a two-stage  furnace, 
the  major  portion  of  the  heating  load  is  supplied  at  the  reduced  input,  with  the  auxiliary  electrical 
components  also  operating  under  the  settings  for  the  reduced  input  rate. 

(3)  For  units  employing  post  purge  after  the  burner  is  shut  off,  the  standard  provides  a 
calculation  procedure  that  limits  the  post  purge  time  to  three  minutes.  For  some  oil-fired  boilers 
that  require  a longer  post  purge  time,  the  formula  in  the  standard  will  produce  erroneous 
(smaller)  flue  loss  result. 

(4)  The  standard  listed  a group  of  integral  functions  named  as  F3  to  F8  that  were  used  to 
produce  the  graphs  in  the  standard  for  hand  calculation  of  the  off-period  losses,  and  in  the 
numerical  algorithm  for  the  losses  in  a computerized  calculation  procedure.  There  is  no 
documentation  on  the  derivation  of  those  functions.  A search  of  the  publications  [3,4]  where  the 
functions  were  first  presented  did  not  show  the  procedure  or  provide  information  about  how  the 
functions  were  derived. 
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This  report  presents  a proposed  revision  to  the  burner  on/off  times  for  a modulating  appliance 
that  corrected  the  inconsistency  problem  in  AFUE,  and  a revised  BOH  calculation  procedure  for 
a modulating  appliance  that  uses  the  measured  fuel  and  auxiliary  electric  inputs  at  the  actual 
operating  conditions.  The  annual  fuel  and  electric  energy  consumption  are  calculated  on  the  basis 
of  the  BOH.  A procedure  for  the  post  purge  calculation  beyond  three  minutes  is  derived,  and  the 
background  information  of  the  process  in  the  derivation  of  the  integral  functions  F3  to  F6  is 
shown. 

2.  ON/OFF  TIMES  FOR  MODULATING  FURNACE/BOILER 

2.1  Background 

The  current  test  procedure  in  ASHRAE  103-1993  prescribes  the  calculation  of  the  AFUE  of  a 
single  stage  furnace  by  using  test  data  from  the  steady  state,  cool  down  and  heat  up  tests  at  the 
rated  input.  In  the  calculation  procedure,  the  average  furnace  cycling  on/off  times  are  assumed 
to  be  3.87  min  on  and  13.3  min  off  based  on  the  characteristic  responses  of  the  room  thermostat 
and  a furnace  oversize  factor  of  0.7  with  respect  to  the  building  heating  load.  This  results  in  the 
furnace  operating  at  an  average  of  22.5  % of  the  rated  capacity  during  the  heating  season.  The 
AFUE  of  a two-stage  furnace  is  calculated  by  using  test  data  from  the  steady  state,  cool  down 
and  heat  up  tests  at  both  the  maximum  (high)  input  rate  and  at  the  minimum  (low)  input  rate. 
There  the  average  furnace  cycling  on/off  time  are  assumed  to  be  10  min  on  and  10  min  off.  And 
this  assumes  the  two-stage  furnace  operating  at  an  average  of  50  % of  the  rated  capacity  during 
the  heating  season. 

During  the  past  years,  a discrepancy  in  the  calculation  procedure  was  discovered  with  respect  to 
some  modulating  furnaces.  The  AFUE  value  for  the  furnaces  operating  as  a two-stage  furnace 
was  lower  than  the  value  when  the  same  test  data  were  used  in  a calculation  with  the  furnace 
operating  as  a single  stage  furnace  operating  under  low  firing  rate.  These  furnaces  were  non- 
condensing two-stage  furnaces  with  a low  fire  to  high  fire  ratio  of  0.60  to  0.70  that  are  typical  of 
two-stage  furnaces.  The  difference  in  the  AFUE  values  between  the  two  calculations 
(modulating  procedure  versus  single  stage  procedure)  can  be  more  than  1 percentage  point. 
However,  based  on  the  test  procedure,  a two-stage  furnace  with  a ratio  of  low  to  high  firing  rates 
of  0.65  would  operate  in  the  low  firing  rate  to  meet  more  than  95  % of  the  annual  heating  load  in 
a heating  season.  Therefore,  the  resulting  AFUE  was  expected  to  be  close  to  the  value  calculated 
on  the  basis  of  single-stage  firing  at  the  low  rate  over  the  entire  heating  season.  The  question  was 
raised  as  to  the  reason  for  this  one  percentage  point  discrepancy  in  the  AFUE  values. 

2.2  AFUE  Calculated  by  Current  Test  Procedure 

Using  test  data  of  a two-stage  furnace,  NIST  investigated  the  cause  of  this  discrepancy  using  the 
computer  program  AFUEBF  (a  computerized  version  of  the  calculation  procedure  of  ASHRAE 
103-1993).  The  original  version  and  a slightly  modified  version  of  the  program  were  run  using 
the  same  test  data  as  input.  The  modified  version  changed  (only)  the  furnace  on/off  times  to 
examine  the  effect  of  the  changes  in  the  on/off  times.  The  changes  made  on  the  on/off  times 
were  (1)  3.87/13.3  min  for  both  single  stage  and  two-stage  calculation,  and  (2)  10/10  min  for 
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both  single  stage  and  two-stage  calculation.  It  was  found  that  the  discrepancy  in  the  AFUE 
described  above  was  caused  by  the  difference  in  the  on/off  times  assigned  to  a single  stage 
furnace  versus  to  a two-stage  furnace.  The  following  examples  show  the  results  (AFUE)  of  two 
gas-fired,  two-stage  non-condensing  furnaces  using  test  data  from  a manufacturer.  Both  have  a 
low  fire  to  high  fire  input  ratio  of  0.65.  The  maximum  input  of  Example  1 is  80  000  Btu/h.  The 
maximum  input  of  Example  2 is  60  000  Btu/h. 


Example  1: 

AFUE  Value 

Using  Single  Stage  Calculation 

Using  Two-Stage  Calculation 

On/Off 

High-Fire 

Low-Fire 

High-Fire 

Low-Fire 

Weighted 

Time* 

Average 

3.87/13.3 

82.64 

81.99 

82.10 

81.85 

81.85 

10.0/10.0 

82.37 

81.59 

81.50 

81.53 

81.52 

*Note:  Current  on/off  times  are  3.87/13.3  min  for  single  stage  furnace,  10/10  min  for  two-stage 
furnace.  Value  in  bold  type  indicates  result  by  current  procedure  for  modulating  furnace. 

Example  2: 

AFUE  Value 

Using  Sing 

;le  Stage  Calculation 

Using  Two-Stage  Calculation 

On/Off 

High-Fire 

Low-Fire 

High-Fire 

Low-Fire 

Weighted 

Time* 

Average 

3.87/13.3 

82.32 

81.31 

81.79 

81.35 

81.35 

10.0/10.0 

81.50 

80.33 

80.34 

80.32 

80.34 

*Note:  Current  on/off  times  are  3.87/13.3  min  for  single  stage  furnace,  10/10  min  for  two-stage 
furnace.  Value  in  bold  type  indicates  result  based  on  current  procedure  for  modulating  furnace. 


Example  1 shows  that  the  current  calculation  procedure  for  a two-stage  furnace  gives  a lower 
AFUE  value  (81.52  %)  than  if  the  furnace  were  operating  as  a single  stage  furnace  (where  the 
on/off  times  are  set  at  3.87/13.3  min)  at  either  high  or  low  firing  rates  (AFUE  of  82.64  % and 
81.99  %,  respectively).  When  the  on/off  times  were  changed  to  10.0/10.0  min  for  the  low-firing, 
single  stage  case,  the  calculated  AFUE  became  81.59  % which  is  closer  to  the  81.52  % for  the 
weighted  two-stage  result  at  the  same  10.0/10.0  on/off  times.  (From  ASHRAE  103-1993,  the 
weighting  factor  for  the  two-stage  calculation  is  greater  than  0.95  at  low  fire  and  less  than  0.05  at 
high  fire  - that  is,  the  two-stage  furnace  operated  almost  all  under  low  fire  to  meet  the  annual 
heating  load).  The  same  close  agreement  happens  if  the  on/off  times  are  3.87/13.3  min  for  both 
single  and  two-stage  calculation  (81.99  % under  single  stage  low  fire  and  81.85  % under  two- 
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stage).  Example  2 is  from  test  data  of  another  (smaller  capacity)  furnace.  The  results  of  this 
example  show  that  the  discrepancy  under  the  current  calculation  procedure  could  be  1 percentage 
point  in  favor  of  the  single  stage  calculation. 

As  described  above,  the  discrepancy  in  the  results  in  the  current  calculation  procedure  is  caused 
by  the  difference  in  the  assigned  on/off  times  for  the  two  types  (single  stage  versus  modulating) 
of  furnaces.  A possible  revision  would  be  to  change  the  assigned  on/off  times  to  some  more 
appropriate  values.  The  problem  is  to  determine  what  should  be  the  corrected  on/off  times  for  a 
modulating  furnace,  or  a modulating  boiler  (where  on/off  time  are  specified  as  15/15  min  in  the 
present  procedure). 

Section  2.3  describes  the  current  procedures  for  determining  on/off  times.  Section  2.4  describes 
the  proposed  procedure  for  determining  the  on/off  times  for  a modulating  fumace/boiler  based 
on  its  minimum  (low-fire)  to  maximum  (high-fire)  capacity  ratio. 


2.3  Review  of  Current  Procedure  for  Single  Stage  Furnaces/Boilers 

In  the  development  of  the  original  test  procedure  for  a single  stage  unit,  the  on  and  off  times  of  a 
fumace/boiler  are  calculated  by  the  following  equations  . 


tON  ~ 


60  X 

4 N X(l-X) 


(2-1) 


60 

t0FF  ~ 4 N X(l-X) 


~ tON 


(2-2) 


where 

N = number  of  burner  cycles  per  hour  based  on  the  responses  of  thermostat  at  half  load, 

= 5 for  furnace 
= 2 for  boiler 

X is  the  load  factor  and  is  defined  as: 


65  - Toa  1 
65  -td  / + « 


(2-3) 


where 

Toa  = the  average  outdoor  temperature.  (Toa  = 42  °F  in  the  test  procedure) 

Td  = the  outdoor  design  temperature.  (To  = 5 °F  in  the  test  procedure) 
a = oversize  factor  = 0.7  for  a single  stage  unit 

For  a single  stage  unit,  X = 0.225,  based  on  the  values  of  Toa,  Td,  and  a in  the  test  procedure. 
This  resulted  in  the  values  of  toisi/toFF  of  3.87/13.3  min  as  given  in  the  test  procedure. 
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2.4  Proposed  Procedure  for  Modulating  Furnaces/Boilers 


For  a modulating  furnace  or  boiler,  the  current  test  procedure  assumed  a load  factor  of  X = 0.5 
which  would  give  the  on/off  times  of  6/6  min  (N  = 5)  for  furnaces  and  15/15  min  (N=2)  for 
boilers.  However,  the  actual  on/off  times  for  furnaces  in  ASHRAE  103-1993  was  set  at  10/10 
instead  of  6/6.  A detailed  description  and  the  rationale  for  the  setting  of  on/off  times  for 
modulating  fumaces/boilers  in  the  current  test  standard  was  presented  in  reference  [5].  However, 
since  a modulating  furnace  (or  boiler)  does  not  always  have  the  specific  minimum  input  rate  that 
would  result  in  X = 0.5  when  the  unit  is  operating  in  the  minimum  input  rate,  the  on/off  times 
actually  vary  depending  on  the  ratio  of  the  minimum  to  maximum  input  rates. 

For  modulating  fumaces/boilers,  ASHRAE  103-1993  defined  a balance  point  temperature  Tc  as 
the  change  over  temperature  for  the  unit  to  switch  from  the  minimum  input  rate  to  the 
modulating  input  or  the  maximum  input  rate.  The  balance  point  temperature  is  defined  as  the 
value  of  the  outdoor  temperature  above  which  the  modulating  equipment  is  capable  of  providing 
the  building  heating  load  requirement  at  the  minimum  or  low  input  rate.  For  outdoor 
temperature  below  the  balance  point  temperature,  the  unit  will  switch  to  its  modulating  mode  of 
operation  for  step-modulating  units  or  to  its  maximum  input  rate  for  a two-stage  unit.  In 
ASHRAE  103-1993,  The  balance  point  temperature  is  calculated  as, 

Tc=  65  — (60)(l  + a)^^  . (2-4) 

Qout 

Where  Qout.  r/Qout  is  the  minimum  to  maximum  capacity  ratio  of  the  modulating  unit,  a = 0.7 
is  the  oversize  factor  at  the  maximum  input  rate  of  the  modulating  unit.  It  is  seen  that  the  value 
of  Tc  depends  on  the  ratio  of  the  minimum  to  maximum  capacity  ratio.  With  a = 0.7  and  a 
capacity  ratio  of  approximately  0.5884,  Tc  would  be  equal  to  approximately  5 °F,  the  average 
outdoor  design  temperature.  This  means  that  for  a fumace/boiler  with  a capacity  ratio  of  greater 
than  or  equal  to  0.5884,  the  fumace/boiler  will  operate  in  the  minimum  (low)  firing  rate  for  the 
entire  heating  season.  For  values  of  the  capacity  ratio  less  than  0.5884,  Tc  would  have  a value  of 
between  5 °F  and  65  °F.  In  this  range,  the  modulating  unit  will  operate  in  the  low  firing  rate 
when  the  outdoor  temperature  is  between  Tc  and  65  °F.  The  unit  will  switch  to  the  maximum 
(high)  firing  rate  (two-stage  control)  or  to  the  intermediate  (step-up)  firing  rates  (step-modulating 
control)  when  the  outdoor  temperature  is  between  the  design  temperature  of  5 °F  and  Tc- 
ASHRAE  103-1993  provides  graphs  to  determine  the  fractions  of  the  annual  heating  load  Xr  and 
Xh  (as  functions  of  Tc)  that  will  be  supplied  by  the  modulating  equipment  operating  at  the 
minimum  (low)  and  the  maximum  (high)  firing  rates,  respectively.  The  standard  also  gives  the 
values  of  the  average  outdoor  temperatures  Toa,  r (between  Tc  and  65  °F  in  low  firing  rate 
operation)  and  Toa.  h (between  5 °F  and  Tc  in  high  firing  rate  operation),  as  function  of  Tc- 
These  temperatures  are  equivalent  to  the  Toa  of  42  °F  for  the  single  stage  equipment. 

For  the  values  of  the  capacity  ratio  (low  fire  to  high  fire)  greater  than  approximately  0.59,  the 
modulating  fumace/boiler  operates  mostly  in  the  low  firing  rate.  If,  as  assumed,  the  oversize 
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factor  a of  the  unit  operating  at  the  high  firing  rate  is  at  0.7,  an  oversize  factor  ocr  at  the 
minimum  firing  rate  with  respect  to  the  building  design  heating  load  can  be  defined  as, 


1 + a* 


Qqut.r 

DHR 


where, 

OCr  = the  oversize  factor  at  the  reduced  input  rate, 

DHR  = Design  Heating  Requirement. 

As  stated  before,  at  the  maximum  firing  rate,  the  oversize  factor  a is  defined  as. 


l + a=&2- 

DHR 

Combining  the  above  two  expressions  results  in, 

7 / 7 \ OUT  .R  f 

l + aR=(l  + a)  — (2-5) 

&OUT 

It  was  noted  at  the  meetings  of  the  multi-fire  working  group  of  ASHRAE  SPC  103R  that  at  the 
values  of  a = 0.7  and  Qout.r  / Qout  of  approximately  0.59,  the  unit  is  no  longer  oversized  at  the 
minimum  firing  rate  (ocr  = 0).  Therefore,  equation  (2-5)  holds  for  values  of  a and  the  capacity 
ratio  that  give  a value  of  (Xr  equal  or  greater  than  0. 

With  the  calculated  ocr  from  equation  (2-5),  equation  (2-3)  for  the  load  factor  A can  be  used  to 
calculate  a load  factor  Xr  at  the  minimum  input  and  a load  factor  Xh  at  the  maximum  input  rate, 
as. 


<55  - Toa.r  1 
65  -Tc  l + 


( ccR  >0.) 


(2-6) 


65  - T oa.h  7 
65  - Td  1 + a 


(2-7) 


It  was  noted  at  the  meetings  of  the  multi-fire  working  group  of  ASHRAE  SPC  103R  that  the 
minimum  outdoor  temperature  for  the  load  factor  Xr  is  Tc,  not  the  design  temperature  Td,  since 
Tc  is  the  limiting  low  outdoor  temperature  for  the  modulating  unit  to  operate  in  low  firing  rate. 

In  ASHRAE  103-1993  only  the  on/off  times  based  on  the  low  firing  rate  is  used  in  the 
calculation  procedure.  This  is  because  that,  for  the  majority  of  modulating  units,  the  requirement 
of  optimum  fuel  to  air  ratio  under  all  firing  rates  results  in  the  capacity  ratio  of  low  to  high  firing 
of  from  0.4  to  0.7.  In  this  range,  from  80  % to  99%  of  the  annual  heating  load  are  met  by  the 


6 


modulating  unit  operating  at  low  fire.  Employing  the  same  approach,  equation  (2-6),  in  place  of 
equation  (2-3),  can  be  used  to  calculate  the  load  factor  X.  The  on/off  times  of  equations  (2-1)  and 
(2-2)  for  the  modulating  unit  can  then  be  calculated  based  on  the  value  of  X from  equation  (2-6). 
The  calculated  on/off  times  are  then  used  in  the  calculating  procedure  for  modulating  furnaces  or 
boilers.  The  following  table  shows  the  calculated  on/off  times  for  a modulating  furnace  (with 
N=5)  and  a modulating  boiler  (N=2),  at  capacity  ratios  of  0.4  to  0.75,  where  0.75  is  the  highest 
ratio  defined  in  ASHRAE  103-1993  for  modulating  units.  The  two  columns  in  the  table  labeled 
"(N=3)"  showed  the  values  of  toN  and  toFF  for  the  case  that  the  value  of  N is  assumed  to  be  3 
instead  of  5 for  modulating  furnaces.  The  value  N=3  resulted  from  the  approximation  of  the 
on/off  of  10/10  for  modulating  furnaces  in  ASHRAE  103-1993: 

Modulated  Furnace  On/Off  Times  Based  on  Quit  r / Pout 


Capacity  Tc 
Ratio  °F 

1+OtR 

Toa.r* 

°F 

Toa.h* 

°F 

XR* 

X=Xr 

Load 

Factor 

tON 

(min) 

(N 

tOFF 

(min) 
= 5) 

ii 

o 2 

4—j  e 

tOFF 

(min) 

3) 

toN 
(min) 
(N  = 

Uf  o 
§■3 

0.40 

24.2 

1.0 

43.5 

19 

0.83 

0.527 

6.34 

5.69 

10.57 

9.49 

15.85 

14.23 

0.45 

19.1 

1.0 

42.6 

15 

0.89 

0.488 

5.86 

6.15 

9.77 

10.25 

14.65 

15.37 

0.50 

14.0 

1.0 

41.8 

11 

0.94 

0.455 

5.50 

6.59 

9.17 

10.99 

13.76 

16.49 

0.55 

8.9 

1.0 

41.0 

7.5 

0.97 

0.428 

5.24 

7.01 

8.74 

11.69 

13.11 

17.53 

0.60 

5 

1.02 

40.7 

5 

0.99 

0.397 

4.98 

7.56 

8.29 

12.59 

12.44 

18.89 

0.65 

5. 

1.11 

40.7 

5 

0.99 

0.367 

4.74 

8.19 

7.89 

13.64 

11.84 

20.46 

0.70 

5 

1.19 

40.7 

5 

0.99 

0.340 

4.55 

8.81 

7.58 

14.69 

11.37 

22.04 

0.75 

5 

1.26 

40.7 

5 

0.99 

0.318 

4.40 

9.44 

7.33 

15.74 

10.99 

23.61 

Single  Stage  (specified  in  ASHRAE 

103) 

0.225 

3.87 

13.33 

— 

— 

9.68 

33.26 

*Note:  The  values  of  Toa,  r (average  outdoor  temperature  at  minimum  input),  Toa,  h (average 
outdoor  temperature  at  maximum  input),  and  XR  (the  fraction  of  annual  heating  load  supplied  at 
the  minimum  input)  as  function  of  Tc  are  from  Figures  23  and  24  of  ASHRAE  103-1993  [2]. 


For  comparison,  the  last  row  in  the  above  table  listed  the  values  of  the  load  factor  X and  the 
burner  on/off  times  for  a single  stage  furnace  (with  N = 5)  and  a single  stage  boiler  (with  N = 2) 
as  specified  in  ASHRAE  103-1993. 

2.5  Examples  of  AFUE  Calculated  on  the  Basis  of  Proposed  Procedure 

The  following  are  two  examples  of  using  the  calculated  on/off  times  (with  N = 5)  in  the  AFUE 
calculation  for  the  two-stage  furnace  examples  shown  previously.  In  place  of  the  fixed  on/off 
values  of  10/10  (minutes)  for  the  furnace,  the  values  were  calculated  using  equations  2-1,  2-2, 
and  2-6)  on  the  basis  of  the  capacity  ratio  Qout,  r /Qout-  The  capacity  ratio  was  computed  based 
on  data  from  the  steady  states  test. 


7 


Example  1: 


AFUE  Value 


Single  Stage  Calculation 

Two-Stage  Calculation 

On/Off 

High-Fire 

Low-Fire 

High-Fire 

Low-Fire 

Weighted 

Times  (min) 

Test  Data 

Test  Data 

Test  Data 

Test  Data 

Average 

3.87/13.3 

82.64 

81.99 

10.0/10.0 

81.50 

81.53 

81.52 

4.59/8.67  (Calculated) 

81.90 

81.70 

81.71 

Example  2: 

AFUE  Value 

Single  Stage  Calculation 

Two-Stage  Calculation 

On/Off 

High-Fire 

Low-Fire 

High-Fire 

Low-Fire 

Weighted 

Times  (min) 

Test  Data 

Test  Data 

Test  Data 

Test  Data 

Average 

3.87/13.3 

82.32 

81.31 

10.0/10.0 

80.34 

80.32 

80.34 

4.58/8.69  (Calculated) 

81.38 

80.91 

80.94 

In  example  1,  the  calculated  AFUE  is  81.99  % based  on  single-stage  at  low  fire  at  the  existing 
on/off  times,  and  81.71  % based  on  two-stage  with  the  new  on/off  times.  For  the  second  furnace, 
the  comparison  is  81.31  % to  80.94  %.  It  is  noted  that  even  though  the  differences  are  smaller, 
the  results  still  do  not  agree  with  each  other  exactly.  This  is  because  that  for  single-stage  low- 
fire,  the  on/off  times  really  should  not  be  3.87/13.3  but  instead  should  be  approximately  4. 6/8. 7 
since  at  the  low  firing  rate,  the  oversize  factor  is  not  0.7,  but  a smaller  value  of  approximately  ccr 
= 0.1 1 (see  previous  table  on  the  value  of  (1+  ocr  = 1.11)  at  Qout.r/  Qout  of  0.65).  If  the  on/off 
values  of  4. 6/8. 7 are  used  in  the  single-stage  calculation,  the  result  would  be  very  close  to  the 
two-stage  procedure  since  in  the  two-stage  operation,  99  % of  the  annual  heating  load  is  met  by 
the  furnace  operating  under  low-fire. 

Lastly,  it  is  noted  that  in  the  above  proposal,  the  on/off  times  based  on  the  low-fire  condition  are 
used  for  the  high-fire  portion  of  the  calculation  also.  This  is  a simplification  based  on  the  reason 
that  the  fraction  of  high  fire  condition  is  small  in  the  practical  range  of  the  capacity  ratios  from 
0.40  to  0.75.  If  more  precision  in  the  result  is  desired,  separate  load  factors  for  Xh  and  Xr  can  be 
calculated  from  equations  (2-6)  and  (2-7),  and  they  in  turn  can  be  used  to  calculate  the  separate 
on/off  times  for  low  fire  and  high  fire  conditions  from  equations  (2-1)  and  (2-2).  However,  as 
stated  before,  it  is  considered  at  the  present  time  to  be  not  necessary  for  the  range  of  capacity 
ratios  encountered  in  practice. 
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3.  ANNUAL  BURNER  OPERATING  HOURS  (BOH)  AND  ANNUAL  ENERGY 

CONSUMPTION  OF  MODULATING  FURNACES  AND  BOILERS 


3.1  Background 

In  the  existing  ASHRAE  103-1993/DOE  test  procedure  [1,2],  an  equivalent  number  of  burner 
operating  hours  (BOHss)  based  on  operation  condition  at  the  maximum  rated  input  was  defined 
to  estimate  the  annual  fuel  energy  consumption  of  both  the  single-stage  and  the  modulating 
fumaces/boilers.  For  single  stage  units,  the  formulation  is  straightforward.  The  number  of  burner 
operating  hours  is  calculated  as  the  annual  building  load  divided  by  the  useful  outputs  of  the 
furnace/boiler  delivered  to  the  conditioned  space.  The  useful  output  energy  from  the 
furnace/boiler  composed  of  the  heat  output  from  the  burning  of  fuel  at  the  heating  seasonal 
efficiency  plus  the  heat  from  the  electrical  power  input  to  the  auxiliary  electric  equipment.  The 
annual  fuel  energy  consumption  is  calculated  as  the  number  of  burner  operating  hour  times  the 
rated  energy  input  plus  the  off-period  pilot  light  consumption  if  employed.  For  modulating 
furnace/boiler,  the  same  formulation  is  retained.  However,  in  the  BOHss  formula,  a 
multiplication  factor  R = 2.3  or  3.0  (depending  on  the  capacity  ratio  of  minimum  (low-fire)  to 
maximum  (high-fire))  is  introduced  in  the  auxiliary  electrical  power  term.  The  heat  energy 
output  from  the  fuel  is  calculated  as  the  maximum  rated  input  at  the  weighted  heating  seasonal 
efficiency  of  the  modulating  furnace/boiler.  The  purpose  of  the  factor  R was  to  account  for  the 
greater  on-time  fraction  per  cycle  (0.50)  for  modulating  unit  as  compared  to  the  on-time  fraction 
per  cycle  (0.225)  for  single  stage  unit.  The  ratio  of  0.50  to  0.225  results  in  R as  2.30.  The 
assumption  was  that  the  auxiliary  electrical  energy  consumption  of  the  modulating  furnace/boiler 
due  to  the  assumed  larger  fraction  on  time  per  cycle  would  be  approximated  by  the 
multiplication  factor  2.3  to  the  auxiliary  electrical  power  input  of  air  circulation  blowers  and 
power  burners,  etc.  Since  the  maximum  fuel  input  rate  and  the  steady-state  auxiliary  electrical 
power  input  measured  at  the  maximum  firing  rate  are  used  to  compute  the  BOHss,  the  term  no 
longer  represents  the  actual  annual  burner  operating  hours  of  the  equipment  which  mostly 
operated  in  the  reduced  or  modulating  input  mode.  Rather,  it  is  used  as  an  equivalent  BOHss 
term  to  approximate  or  estimate  the  annual  fuel  energy  consumption  of  the  modulating 
equipment  during  a heating  season.  (Appendix  C of  ASHRAE  103-1993.)  However,  these 
assumptions  ignored  the  effects  of  the  possible  different  values  of  the  electrical  power  input  at 
the  maximum  and  reduced  fuel  input  rates  for  furnaces/boilers  with  a variable  speed  circulating 
air  blower  and  combustion  blower.  Also,  the  use  of  the  maximum  fuel  input  rate  in  the 
calculation  formula  for  the  BOHss  causes  confusion  to  users  of  the  standard  since  for  a 
modulating  furnace/boiler,  the  design  purpose  is  for  the  equipment  to  operate  at  the  reduced 
input  rate  most  of  the  heating  season. 

To  improve  upon  the  above-mentioned  assumptions  in  the  calculation  procedure,  the  proposed 
calculation  procedure  would  use  the  actual  power  inputs  of  a variable  speed  blower  at  the  firing 
rates  it  is  operating  under.  The  actual  burner  operating  hours  at  the  maximum  input  and  the 
reduced  (minimum)  input  rates  are  calculated  to  determine  the  annual  energy  consumption  of  the 
equipment  operating  under  both  modes.  The  procedure  would  also  use  both  the  maximum  and 
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the  reduced  fuel  input  rates  in  the  calculation  of  the  burner  operating  hours  under  each  firing  rate 
to  eliminate  the  approximation  of  an  equivalent  BOHss- 

3.2  Review  of  Test  Procedure  for  Single  Stage  Furnaces/Boilers 

From  Appendix  A of  NBSIR  78-1543  [3]  or  Appendix  C of  ASHRAE  103-1993  [2],  the  number  of 
Burner  Operating  Hours  (BOH)  and  the  annual  fuel  consumption  Ep  for  a single  stage  furnace 
are  calculated  as,  for  indoor  installation, 

BOH  [Heating  Season  Fuel  Consumed  - (Qp)(Heating  Season  Hour)]  / Qin, 

Ef  =Annual  Fuel  Energy  Consumption 

= Heating  Season  Fuel  Consumed  +,(Qp)(Non-Heating  Season  Hours), 

Heating  Season  Fuel  Consumed  = [Annual  Building  Load  - (PE  + y BE)(BOH)]  / Effyns, 

Annual  Building  Load  = (C)(HLH)(DHR, 

where, 

Qin  = steady-state  rated  fuel  energy  input 
Qp  = pilot  light  fuel  consumption  rate 
BE  = electric  power  to  circulation  blower  or  water  pump 
PE  = electric  power  to  burner 

y = ratio  of  circulation  blower  or  pump  on-time  to  burner  on-time  in  a burner  cycle 

EffyHs  = heating  seasonal  efficiency 

HLH  = Heating  Load  Hours 

DHR  = design  heating  requirement  of  a building 

C = 0.77,  an  adjustment  factor  defined  in  Appendix  C of  ASHRAE  103-1993. 

For  uniform  rating  purpose,  ASHRAE  103-1993  established  the  following  national  average 
values  for  some  of  the  weather-related  parameters: 

Non-Heating  Season  Hours  = 4600, 

Heating  Season  Hours  = 4160, 

Heating  Load  Hours  (HLH)  = 2080. 

For  simplicity,  let  Ep  = [(PE)  + (y)(BE)],  rju  = Effyns,  and  using  the  national  average  values, 

BOH  = HLH  [(C)(A)(DHR)  - B]  = 2080  [(C)(A)(DHR)  - B],  (3-1) 

A = 1 / [Ep  + (Qin  - QP)  T|u], 

B = (Qp)(r|u)(A)(Heating  Season  Hours)  / HLH  = 2 Qp  r|u  A. 

The  above  equations  are  the  same  as  those  given  in  Appendix  C of  ASHRAE  103-1993,  except 
that  for  simplicity,  the  percent  and  unit  conversion  factors  (100  in  percent  and  3413  W = 1 Btu/h) 
spelled  out  in  the  equations  in  ASHRAE  103-1993  are  not  included  here.  (For  use  in  actual 
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calculation  with  numbers,  consistency  in  units  with  the  appropriate  unit  conversion  factors  is 
needed  and  included  in  the  equations  in  the  standard.) 

3.3  Proposed  Procedure  for  Modulating  Furnaces/Boilers 

3.3.1  Annual  Burner  Operating  Hours  (BOH): 

The  following  nomenclature  is  used  in  the  derivation  of  the  proposed  procedure: 

BQHh  = annual  BOH  of  a two-stage  unit  at  the  maximum  rated  input  Qin.h 
Qin.h  = fuel  input  rate  at  the  maximum  firing  rate 

EpH  - electric  power  consumption  measured  under  steady  state  at  input  Qin.h 
rjUiH  = part-load  efficiency  of  a two-stage  unit  operating  at  input  Qin.h 

BOHm  = annual  burner  operating  hours  of  a step-modulating  unit  at  modulating  mode 
Qin.m  = fuel  input  rate  in  the  modulating  mode  of  a step-modulating  unit  as  calculated  in  sec.  1, 
Appendix  C of  ASHRAE  103-1993 

EpM  - electric  power  input  of  a step-modulating  unit  operating  in  the  modulating  mode  and 
measured  under  steady  state  at  input  Qin.h 

riu.M  = average  part-load  efficiency  of  a step-modulating  furnace/boiler  operating  in  the 
modulating  mode  (Sec.  11.4.9.2.3,  ASHRAE  103-1993) 

BOHr  = annual  burner  operating  hours  at  the  minimum  input  Qin.r 
Qin.r  = fuel  input  rate  at  the  minimum  firing  rate 
EpR  = electric  power  input  measured  under  steady  state  at  input  Qin.r 
rju.R  = part-load  efficiency  operating  at  input  Qin.r 

It  is  assumed  that  for  the  modulating  furnace/boiler  operating  at  either  the  maximum  or 
minimum  input,  or  at  modulating  mode,  the  Design  Heating  Requirement  (DHR)  of  the  building 
is  not  changed  and  therefore  the  required  fumace/boiler  heat  output  to  meet  the  DHR  of  the 

building  is  not  changed. 

For  simplicity  in  showing  the  derivation  process,  it  is  assumed  that  the  modulating  unit  has  an 
electric  ignition  device  so  that  the  pilot  Qp  = 0 (extension  of  results  to  include  Qp  will  be  shown 
later  in  this  report).  With  this  assumption,  the  general  equations  for  the  BOH  (Eq.3-1)  and  the 
variables  A and  B of  Section  3.2  above  are  reduced  to: 

B = 0, 

A = 1 / [Ep  + (Qin)  rju] 
and 

BOH  - 2080  [(C)(A)(DHR)].  (3-2) 

With  the  above  assumptions,  the  burner  operation  hours  in  meeting  the  total  annual  heating  load 
under  the  maximum  firing  rate  is. 
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BOHh*  = HLH  [(C)(Ah)(DHR)], 

Ah  = 1 / [Epn  + (Qin.h)  Tu.h] 

where  ihe  asterisk  * indicates  the  total  BOH  if  the  appliance  operates  at  a particular  (maximum 
here)  firing  rate  for  the  entire  heating  season.  If  only  a fraction  of  the  total  annual  heating  load  is 
met  by  the  appliance  operating  at  the  maximum  firing  rate,  the  actual  BOH  at  the  maximum 
firing  rate  (for  a two-stage  unit)  during  the  heating  season  will  be, 

BOHh  = (XH)  (BOHh  ) = (Xh)[(C)(HLH)(Ah)(DHR)],  (3-3) 

where  XH  is  defined  in  Sec.  1 1.4. 8. 6 of  ASHRAE  103-1993  as  the  fraction  of  the  total  annual 
heating  load  that  is  provided  by  a two-stage  appliance  at  the  maximum  firing  rate. 

Similarly,  the  total  burner  operating  hours  to  satisfy  the  total  annual  heating  load  when  operating 
under  the  minimum  firing  rate  for  the  entire  heating  season  is, 

BOHr*  = HLH  [(C)(Ar)(DHR)], 

Ar  = 1 / [EpR  + (Qin.r)  T|u.r]. 

And,  if  only  a fraction  of  the  total  annual  heating  load  is  met  by  the  appliance  operating  at  the 
minimum  firing  rate,  the  BOH  at  the  minimum  firing  rate  during  the  heating  season  will  be, 

BOHr  = (XR)  (BOHr*)  = (Xr)[2080  (C)(Ar)(DHR)],  (3-4) 

where  Xr  is  defined  in  Sec.  11.4.8.7  of  ASHRAE  103-1993  as  the  fraction  of  the  total  annual 
heating  load  that  is  provided  by  a modulating  appliance  at  the  minimum  firing  rate. 

For  a step-modulating  furnace/boiler,  the  fuel  input  rate  at  the  modulating  mode  is  not  measured 
in  ASHRAE  103-1993  since  the  firing  rate  at  the  modulating  mode  varies  and  cannot  be 
determined  a priori.  Instead,  ASHRAE  103-1993  calculates  a Qcm.m  as, 

Qin.m  = Qout.m  / Effyss.M,  (From  Sec.  1 of  Appendix  C,  ASHRAE  103-1993) 
where 

Effyss,M=  average  steady-state  efficiency  calculated  in  Sec.  1 1.4. 8. 8 of  ASHRAE  103-1993 
Qout.m  = average  heat  output  rate  calculated  in  Sec.  1 1.4.8.10  of  ASHRAE  103-1993 

The  fraction  of  the  total  annual  heating  load  that  is  provided  by  a step-modulating  appliance  in 
the  modulating  mode,  XM,  is  assumed  in  ASHRAE  103-1993  to  equal  (1-  Xr).  That  is,  the 
modulating  mode  is  at  any  input  rate  greater  than  the  minimum  input  rate.  Also,  the  electrical 
power  inputs  to  the  step-modulating  appliance  in  the  modulating  mode  (BEm  and  PEm)  are  not 
measured  in  ASHRAE  103-1993  and  are  assumed  to  be  those  measured  during  the  steady-state 
test  at  the  maximum  fuel  input  rate,  that  is,  EpM  = Epn-  Therefore,  the  burner  operating  hours  for 
a step-modulating  appliance  operating  in  the  modulating  mode  during  the  heating  season  is, 


12 


BOHm  - (Xh)[2080  (C)(Am)(DHR)], 


(3-5) 


where 

Am  = 1 / [Epn  + (Qin.m)  t1u,m] 

In  the  above  formulation,  the  subscripts  H,  M and  R denote  the  values  at  the  maximum  firing 
rate,  the  modulating  mode  (step-modulating),  and  the  minimum  firing  rate,  respectively.  The 
variables  Epn,  EpR,  Qin,h>  Qin.r,  Qin.m,  “Hu.h,  Tu.r  and  tVm  are  defined  previously. 


3.3.2  Annual  Fuel  Energy  Consumption 

For  two-stage  modulating  appliances,  the  annual  heating  load  is  met  by  the  appliance  operating 
at  either  the  minimum  firing  rate  or  at  the  maximum  firing  rate.  Therefore,  the  annual  fuel  energy 
consumption  for  a two-stage  appliance  is, 

Ef  = (BOHhXQin.h)  + (BOHr)(Qin,r).  (3-6) 

For  step-modulating  appliance,  the  total  annual  heating  load  is  met  by  the  appliance  operating  at 
either  the  minimum  firing  rate  or  at  the  modulating  mode  (up  to  the  maximum  firing  rate). 
Therefore,  the  annual  fuel  energy  consumption  for  a step-modulating  appliance  is, 


Ef  = (BOHmXQin.m)  + (BOHr)(Qin,r). 


(3-7) 


3.3.3  Annual  Auxiliary  Electrical  Energy  Consumption 

The  annual  auxiliary  electrical  energy  consumption  for  a two-stage  appliance  is, 

Eae  = (BOHnXEpn)  + (BOHr)(EPr).  (3-8) 

The  annual  auxiliary  electrical  energy  consumption  for  a step-modulating  appliance  is, 

Eae  = (BOHmXEph)  + (BOHr)(Epr).  (3-9) 

In  the  general  form  without  the  subscript,  Ep  = (BE)  + (y)(PE),  as  defined  previously. 

3.3.4  Units  with  Standing  Pilot  Light 

It  is  assumed  previously  that  the  modulating  appliance  employs  interrupted  electronic  ignition 
and  no  pilot  light  is. needed.  If  the  modulating  appliance  has  a continuing  burning  pilot  light  with 
fuel  input  rate  Qp,  the  above  formulation  should  be  modified  by  setting  the  parameters  Ah,  Ar, 
Am  and  BH,  BR,  BM  as 


Ah  - 1 / [Epn  + (Qin.h  - Qp)  t1u.h], 

Ar  = 1 / [EPr  + (Qin.r  - Qp)  T)u,r], 
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Am  = 1 / [EpH  + (Qin.m  - Qp)  Pu,m] 

Bh  = - Qp  Hu.h  Ah, 

Br  = 2 QP  r|u,R  Ar, 

Bm  = 2 Qp  rju,M  Am 

And  the  burner  operation  hours  will  be  modified  as 

BOHh  = (XhXHLH)  [(C)(Ah)(DHR)  - BH],  (3-10) 

BOHr  = (Xr)(HLH)  [(C)(Ar)(DHR)  - Br],  (3-11) 

BOHm  = (XhXHLH)  [(C)(Am)(DHR)  - BM]  (3-12) 

The  annual  fuel  energy  consumption  will  be  modified  as, 

For  two-stage  modulating  appliance, 

Ef  = (BOHhXQin.h)  + (BOHr)(Q,n.r)  + [8760  - (BOHh  + BOHr)]  Qp  (3-13) 

For  step-modulating  appliance, 


Ef  - (BOHmXQin.m)  + (BOHrXQin.r)  + [8760  - (BOHm  + BOHr)]  Qp  (3-14) 

It  is  assumed  that  the  value  of  Qp  is  the  same  for  all  modes  of  operations.  The  general  equation 
forms  of  the  parameters  A and  B and  the  BOH  are  shown  in  section  3.2  previously. 

Note  that  the  above  formulation  is  for  furnace/boiler  installed  indoors.  For  furnace/boiler 
installed  in  an  unconditioned  space  as  an  isolated  combustion  installation  or  installed  outdoors, 
the  current  DOE  test  procedure  specified  the  procedure  for  modifying  the  calculation  of  the 
burner  operating  hours.  The  modification  involved  the  accounting  of  the  portions  of  the  electric 
energy  input  to  the  blower  motors  that  contribute  to  the  heating  of  the  conditioned  space  and  the 
portions  that  are  considered  as  lost  to  the  ambient. 

4.  CALCULATION  OF  OFF-PERIOD  LOSSES  FOR  OIL-FIRED  BOILERS  WITH 

LONG  POST  PURGE  TIME 


4.1  Background 

As  discussed  previously,  the  post  purge  period  in  the  test  and  calculation  procedures  of 
ASHRAE  103-1993  was  limited  to  three  minutes  after  the  burner  is  shut  off.  This  limit  was 
considered  adequate  for  most  fossil-fueled  appliances.  However,  for  certain  oil-fired  boilers  the 
operating  requirement  is  such  that  the  post  purge  period  runs  longer  than  three  minutes.  Since  the 
longer  post  purge  period  resulted  in  greater  off-period  losses,  it  was  deemed  appropriate  that 
calculation  for  the  longer  post  purge  period  should  be  specified  to  account  for  the  greater  losses. 
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4.2  Calculation  Procedure  for  Losses  during  Post  Purge  Period 


4.2.1  Assumptions 

The  following  assumptions  were  made  during  the  derivation  of  the  post  purge  period  losses  in 
ASHRAE  103-1993.  They  are  used  in  the  present  derivation  also. 

(a)  Constant  volumetric  flue  gas  flow  rate  during  post  purge  to  equal  to  the  on-cycle  rate. 

(b)  Flue  gas  temperature  approximated  by  piecewise  linear  variation  over  small  time  interval  At. 

(c)  Measurement  data  of  the  flue  gas  temperature  in  the  off-period,  TF,oFF(t)  at  time  interval  At 
following  burner  shut-off. 

4.2.2  Calculation  for  Sensible  Loss 

(a)  Assuming  constant  combustion  blower  or  draft  inducer  volumetric  flow  rate  during  the  post 
purge  period  to  equal  the  steady  state  on-penod  rate  results  in: 

^ F.ON  = M F .ON  ! P F.SS  = ^ F .OFF  ^ P F .OFF  = ^F.OFF  (4  ~ 1) 


where  VFi0n  and  Mf,on  are  the  volumetric  and  mass  flow  rate  of  the  flue  gas  during  the  on- 
period,  respectively.  Vf,off  is  the  off-period  volumetric  flow  rate  during  post  purge,  and  pF,ss  and 
Pf.off  are  the  density  of  flue  gas  at  temperatures  TF,ss  (steady  state)  and  TF,oFF(t),  respectively. 
MF,oFF(t)  and  TF,oFF(t)  are  the  flue  gas  mass  flow  rate  and  temperature  at  time  t during  the  post 
purge  period,  respectively. 

From  the  above  equation,  and  assuming  the  flue  gas  as  a perfect  gas. 


M f. OFF  (?) 


M f, ON  P F,  OFF 
P F.SS 


M F,0N  (Tf  ss  + 460) 

(Tf.off  + 46°i 


(4-2) 


(b)  Assuming  piecewise  linear  flue  gas  temperature  variation  over  time  interval  At,  (tj  < t < ti+1), 
during  the  post-purge  period,  gives 

Tf.off  (?)  = (T F 0FF  iti+x  ~TF  0FF  i lti)(—)  — (TF  0FF  i ~TF  0FF  M )(— )(0 

At  At 

where  Tf.off.  i = TF>0FF(h)  = flue  gas  temperature  at  time  tj  in  post  purge  period. 

Applying  the  factor  Ct.off  for  cycling  effect  to  (TF,0FF(t)  - TRA)  as  in  ASHRAE  103-1993: 


Tf.off  (?)  -Tra  ~ ^Y 


OFF 


(Tf.off.  Jm  Tf  off  Mt{ )(  ) T RA 

At 


(Tf.off.i  TF0FFi+ , )(  )(t) 


At' 
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A off  7)  Am  + C, 


OFF 


^f.oFF.i^  + 1 TFQff.  i+lfX  . ) ^ RA 

At 


( Tf.OFF.i  TF'OFF,i+\  X^XOj 

(4-3) 

where  Ct.off  = the  cool-down  temperature  profile  correction  factor  for  the  effect  of  cycling  [1,21. 
TrA  = laboratory  room  air  temperature. 


The  total  sensible  heat  loss  during  the  post-purge  period  is. 


As,  purge  = C P F, OFF  17 F .OFF  A)  ~ Am  J ^ t 


- CPM  F ON  (TFSS  + 460)  J 


f ^ F .OFF  7)  Am 


0 A\ off  A) + 460 


dt 


— C pM  F 0N(TF  ss  +460)  ^ 


1=1 


'•*}  T (t)-T 

f 1F.OFF\l) 


I 


RA 


T 

t,  1 F.OFF 


(0  + 460 


( from  equation  (4  - 2)  for  M F 0FF ) 
dt  (4-4) 


where  tp  = total  length  of  post  purge  time  interval. 

N = tp  / At  = number  of  intervals  of  length  At  within  tp 
Cp  = specific  heat  capacity  of  flue  gas 

The  percent  sensible  loss  (as  percent  of  total  input  during  the  on-cycle)  during  the  post  purge 
period  is,  by  substituting  the  expressions  for  (Tf.offO)  - TrA)  and  Tf.offO)  above  into  Ls,  purge 
and  integrate: 


Ls, OFF.  purge  ^ ^0 


L 


S,  purge 


= 100 


( Qin  ) ?ON 
CpM  F i0N  (Tf  ss  +460)A t ^ 

( QlN  ) ^ ON  1=1 


1 — 


Am +460 


CTS  (LFi0FF  l Ttf _ 0FF J+l 


M) 


■In 


L ra  + 460  + CTS  (Tf  off  j Am  ) 


Am  + 460+C„(7>t 


OFF , i+l 


Am  ) 


(4-5) 

where  Cjs  = Ct.off  for  indoor  combustion  air, 

= Cs'  Ct.off  for  outdoor  combustion  air, 

Cs'  = off-period  correction  factor  for  outdoor  unit  as  defined  in  ASHRAE  103-1993, 
toN  = burner  on-time  per  cycle, 

Qin  = fuel  input  rate. 

100  = conversion  factor  from  decimal  to  percent 


4.2.3  Calculation  for  Infiltration  Loss 


The  total  infiltration  heat  loss,  L|.purge,  during  post  purge  period  is,  by  the  same  way. 
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L, > l>uri!e  -C P^M F0FF  [Tref  Toa  ] dt 

o 

= CpMF'ON  (Tf  ss  + 460)(7^  - Toa  ) X 


1=1 


'i+l  i 

f l— 

] Tfoff  (0  + 460 


(4-6) 


where  Tref  = a reference  temperature  defined  to  be  70  °F  in  ASHRAE  103-1993, 

Toa  = heating  season  average  outdoor  temperature  as  defined  in  ASHRAE  103-1993. 


The  percent  infiltration  loss  (as  percent  of  total  input  during  the  on-cycle)  during  the  post  purge 
period  is.  by  substituting  the  expression  for  Tf.off(U  in  equation  (4-3)  into  L|,  purge,  and  integrate: 
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OFF , purge 


= 100 
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purge 


100 


(Q IN  ) ^ ON 

CPMF'ON(TF_ss+460)At  " 


T -T 

1 REF  1 OA 


(QlN  ) ^ ON 


I .OFF  ( Tp.OFF.i  Tf 
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OFF.  i+l 


i=l 


Tra  + 460  + CT  0FF  (TF  0FF  , ) 

Tra  + 460  + CT  0FF  (TF  OFF  i+1  — ) 

(4-7) 


It  is  noted  that  when  N = 1 ; At  = tp,  and  the  expressions  for  the  sensible  and  infiltration  losses 
during  post  purge  reduce  to  the  ones  in  ASHRAE  103-1993. 

5.  BACKGROUND  OF  THE  INTEGRAL  FUNCTIONS  FOR  CALCULATING  THE 

OFF-PERIOD  LOSSES  IN  ASHRAE  STANDARD  103-1993 

5.1  Background 

In  ASHRAE  103-1993,  the  off-cycle  sensible  and  infiltration  losses  (see  sections  11.2.10.6  and 
1 1.2.10.8  of  the  Standard)  are  calculated  through  the  pair  of  functions  (F3,  F4);  (F5,  F6);  or  (F7, 
F8)  depending  on  fumace/boiler  location  (indoor,  outdoor,  or  isolated  combustion  system)  and 
stack  damper  installation  (with  or  without).  The  functions  are  presented  in  graphical  forms 
(Figures  13  through  19  of  the  Standard)  for  manual  calculation  with  an  electronic  calculator  as 
well  as  in  formula  forms  (in  terms  of  definite  integrals)  presented  in  Appendix  E of  the  standard 
In  the  formulas  for  the  functions  F4,  F6  and  F8,  there  is  a number  100  appearing  as  a divisor  as 
well  as  an  addition  in  the  numerator  and  denominator  of  the  first  integral  of  each  function  (see 
Appendix  E of  ASHRAE  103-1993  for  the  exact  locations).  The  same  100  also  appears  in  the 
formula  for  the  stack  damper  effectiveness  factor  Do  and  in  the  defining  equations  for  the  terms 
G3  and  G4  in  the  formula  for  Do  (see  section  1 1.2.9.18  of  ASHRAE  103-1993).  In  ASHRAE 
103-1993  the  term  100  was  not  defined.  This  brief  report  shows  the  origin  of  the  term  and  the 
derivations  of  the  integral  functions. 

5.2  Derivation  of  the  Integral  Functions 
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For  illustration  purpose,  the  following  derivations  use  the  case  of  a furnace/boiler  using  indoor 
air  for  combustion  and  installed  indoors  without  a stack  damper.  The  applicable  function  pair  is 
F3  and  F4. 


L 


(100)(0.24)(Z)/,  )[1  + (RT  F )(A/F)]  (Tf  ss -Tm+  530) 


1 19  t0FF 


1.56 


S.OFF 


(HHVa  )(t0N ) 


(Tf,SS  7 RA  ) 


0.56 


dt 


IS  1 OFF 

_ ^ S.OFF  J 
? ON  0 


(?F,OFF  Tra  ) 


1.56 


dt 


0 {Tf<off-Tra+  530)119 
(5-1) 


o (Tf.off-Tra+5 30) 1,9 
In  Appendix  A of  the  report  NBSIR  78-1543  [3],  the  off-cycle  sensible  heat  loss  is  given  as: 


where 

ys F.o  — V F.O.X  C , OFF  = C,  OFF  ^ F .OFF  (*3  ) ~ ^ F .OFF  (°°)1  e 3 

Vf.-  = TF.OFF^)~TRA^t  — oo  denotes  reaching  equilibrium  condition  during  off  - period ) 
T f.off  (0  = V f. o ^ + 7^ 

t3  = a specific  time  in  the  off-period  at  which  the  flue  gas  temperature  is  measured, 

Toff  = off-period  time  constant  of  the  unit. 


With  the  above  defining  equations  for  vj/p.o  and  \|/f,°=,'  Ls.off  can  be  written  as, 


JS  ‘ OFF 

^ S.OFF  f 


K l°FF 

^ S.OFF 


^ S.OFF 


ton  0 (VVo*"/Toff +VV-+530) 


1.19 


dt  = -^£5-  \F{t,VF,0,VF.„)dl  (5-2) 


■ON  0 


It  is  seen  that  in  the  above  expression  for  Ls.off  is  a function  of  t with  the  \|/p,o  and  \\f?„  as 
parameters  whose  values  are  determined  from  test  data  for  a particular  furnace/boiler. 

During  the  1970s  when  the  furnace  test  procedure  was  developed,  the  accurate  numerical 
evaluation  of  the  integral  in  the  above  expression  was  a major  problem  without  using  a central 
computer.  Therefore,  it  was  decided  that  graphs  would  be  developed  to  facilitate  the  calculation 
with  the  use  of  an  electronic  calculator.  However,  since  both  \|/f,o  and  \|/f,<»  are  related  to  a 
specific  furnace/boiler  and  will  have  different  values  for  different  furnace/boiler,  and  a three- 
dimensional  graph  with  both  as  parameters  on  a graph  was  not  practical,  a large  number  of 
graphs  would  be  required.  That  is,  a number  of  graphs  with  vj/p.o  as  the  parametric  variable, 
similar  to  the  two  graphs  for  the  functions  F3  and  F4  in  ASHRAE  103-1993,  would  be  needed 
for  each  value  of  \j/f,o°.  One  way  to  reduce  the  number  of  graphs  required  was  to  take  the 
parameter  ij/p.oo  out  of  the  integral.  During  the  evaluation  process,  it  was  found  that  the  function 
in  the  integral,  F(t,  \]/F.o,  \|/f~),  can  be  approximated  as  a linear  function  of  ij/f.oo  for  constant 
values  of  t between  0 and  toFF  (13.3  min  for  furnaces  and  33.3  min  for  boilers  as  assumed  in  the 
test  procedure)  and  for  constant  \|/F,o-  That  is,  it  can  be  assumed  that 

F(t,\f/F  0,y/F'J  = Cl  y/>,„  + C2  (5-3) 
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with  the  conditions  F = F(t,  vj/F.o,  0)  for  \|/F,~  = 0;  and  F = F(t,  \j/Flo,  100)  for  \j/F„  = 100,  where  the 
point  \j/F,oo  =100  was  arbitrarily  chosen  as  one  possible  value  of  off-cycle  flue  gas  temperature 
(differential)  at  equilibrium.  Usually,  \|/F„  = TpoffI00)  - TRA  would  be  much  smaller  than  100. 
Figure  5- 1 shows  a plot  of  the  linear  variation  of  the  function  F as  a function  of  \|/F,TO  between  the 
values  of  0 to  100  for  ranges  of  t = 0 min  to  30  mm  after  the  burner  was  off  for  a boiler  operating 
indoors  with  the  value  of  \|/F,o  determined  by  test  data  based  on  ASHRAE  103-1993.  It  is  seen 
that  the  linear  approximation  is  very  good. 
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Variation  of  the  function  F versus  \j/F  ^ 


With  the  above  boundary  conditions  (F  = F(t,  \|/p.o,  0)  for  i|/F„  = 0;  and  F = F(t,  \|/F,o,  100)  for 
\|/F,oo  = 100),  the  functional  form  of  C|  and  C2  in  the  linear  approximation  can  be  determined  as. 


C2  = F(t,  1 j/f.o,  0) 

Ci  - [F(t,  if/p.o,  100)  - F(t,  1 f/p.o,  0)]  / 100 


Or,  explicitly, 


C2(t,VF,o)  = F(t,\i/FO,0)  = 
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Substituting  the  expressions  of  Q and  C2  into  the  function  F inside  the  integral  for  the  loss  Ls.off 
and  changing  the  variable  of  integration  from  t to  t/ioFF,  result  in  the  expression  for  the 
approximation  of  Ls.off  by  the  functions  F3  and  F4  as  shown  in  section  1 1.2.10.6  of  ASHRAE 
103-1993. 
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It  is  noted  that  since  the  parameter  xj/p  - is  a constant,  it  is  taken  out  of  the  integral  involving  F4 
in  the  expression  for  Ls.off- 

The  function  F5,  F6,  F7,  and  F8  in  ASHRAE  103-1993  were  derived  the  same  way.  It  is  noted 
here  that  in  ASHRAE  103-1993,  the  derivation  of  the  damper  effectiveness  factor  Do  and  the 
defining  equations  for  the  terms  G3  and  G4  in  the  formula  for  Do  (see  section  1 1.2.9.18  of 
ASHRAE  103-1993)  used  the  same  approximation  for  the  off-cycle  losses  (and  the  same  term 
TOO').  Details  of  that  derivation  is  presented  in  the  ASHRAE  Transaction  paper  titled  "The 
Effect  of  a Time  Delayed  Stack  Damper  on  Off-Cycle  Heat  Losses  for  Residential  Heating 
Equipment"  by  C.  Park,  et  al.[61 

As  mentioned  before,  the  value  100  (as  one  possible  temperature  difference  between  the  flue  gas 
temperature  and  room  temperature  at  equilibrium  during  cool  down)  was  selected  arbitrarily  to 
determine  the  slope  of  the  straight  line  in  the  linear  approximation  of  the  function  F(t,  vj/p.o,  Vf,~). 
For  furnaces  the  actual  temperature  difference  would  be  fairly  small  for  unit  with  pilot  light  (and 
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0 for  unit  with  interrupted  electronic  ignition  as  specified  in  the  test  procedure),  and  for  boilers 
the  actual  temperature  difference  would  be  generally  much  less  than  100. 

6.  SUMMARY 

This  report  documented  the  developments  of  proposed  revisions  to  the  ASHRAE  standard  103- 
1993  on  the  method  of  testing  for  residential  central  furnaces  and  boilers.  The  changes  include 
the  following: 

(1)  A revised  procedure  to  calculate  the  average  burner  on/off  times  of  modulating  furnaces  and 
boilers  during  the  cyclic  operation  of  the  burner.  The  procedure  used  the  same  cycles  per  hour  at 
half  load  as  were  assumed  for  a single  stage  unit  (5  cycles  per  hour  for  furnaces  and  2 cycles  per 
hour  for  boilers),  thus  gives  a more  consistent  comparison  with  a single  stage  unit. 

(2)  A revised  method  for  the  calculation  of  the  annual  burner  operating  hours  using  the  measured 
electric  power  input  at  the  actual  fuel  input  rate  to  determine  the  annual  burner  operating  hours 
and  the  annual  fuel  and  electric  energy  consumption  at  the  actual  burner  firing  rates  of  a 
modulating  unit. 

(3)  An  extension  of  the  current  procedure  to  calculate  the  off-period  flue  loss  of  oil-fired 
equipment  with  a longer  post  purge  time.  The  procedure  will  more  accurately  determine  the  off- 
period  flue  loss  that  is  higher  when  the  boiler  is  operating  with  a longer  than  3 min  of  post  purge 
time. 

This  report  also  provides  the  background  information  on  the  derivation  of  the  integral  functions 
listed  in  Appendix  E of  ASHRAE  standard  103-1993  that  were  not  being  formally  documented. 
Those  integral  functions  were  used  for  producing  the  graphs  in  the  standard  for  the  hand 
calculation  of  the  off-period  losses  with  a calculator,  and  also  in  the  numerical  algorithm  in  the 
computerized  calculation  procedure. 
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